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We present a calculation of the production of neutrinos during propagation of ultra-high energy cosmic rays as may 
be produced in astrophysical sources. Photoproduction interactions are modeled with the event generator SOPHIA 
that represents very well the experimentally measured particle production cross sections at accelerator energies. We 
give the fluxes expected from different assumptions on cosmic ray source distributions, cosmic ray injection spectra, 
cosmological evolution of the sources and different cosmologies, and compare them to the Waxman-Bahcall limit on 
source neutrinos. We estimate rates for detection of neutrino induced showers in a km 3 water detector. The ratio of 
the local high energy neutrino flux to the ultra-high energy cosmic ray flux is a crucial parameter in distinguishing 
between astrophysical and cosmological (top-down) scenarios of the ultra-high energy cosmic ray origin. 
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I. INTRODUCTION 



The highest energy cosmic rays are energetic enough 
to have photoproduction interactions on the microwave 
background. These collisions cause energy loss affecting 
the cosmic ray spectrum [jjj - the GZK cutoff. In most as- 
trophysical environments all secondary mesons produced 
in photoproduction interactions decay into 7-rays and 
neutrinos. Shortly after the original papers on the GZK 
cutoff it was suggested [|J that guaranteed fluxes of ultra- 
high energy (UHE) neutrinos will be produced by the 
propagation of UHE cosmic ray (UHECR) protons in the 
Universe. This suggestion was followed by more sophisti- 
cated estimates |3|-[7| that attempted to predict more re- 
alistically the expected neutrino fluxes and relate the de- 
tection of such fluxes to the neutrino cross section at very 
high energy and the then unknown mass of the W boson. 
Hill & Schramm ||||] introduced the cosmological evolu- 
tion of cosmic ray sources and used the measurements 
of the cosmic ray spectrum by the Haverah Park JlQ ] 
and the Fly's Eye [ |ll| experiments to determine mini- 
mum and maximum allowed normalizations for the flux 
of such 'propagation' neutrinos and calculated detection 
rates for different types of detectors. More recent esti- 
mates include the work of Stecker and collaborators [12| , 



Yoshida & Teshima jl3], and Protheroe & Johnson 14 1. 

Meanwhile the world statistics of ultra-high energy cos- 
mic rays has significantly increased p| and, most impor- 
tantly, two events of energy substantially above 10 20 eV 
were detected by the Fly's Eye [H and AGASA Q ex- 
periments. These events suggest that the maximum en- 
ergy of cosmic ray acceleration E max may be significantly 
higher than the previous nominal estimate of 10 20 eV, if 
these events are not a result of the decay of extremely 
massive exotic particles or other exotic processes [[l9| . 

We assume that UHECR are of astrophysical origin 
and present here a new estimate of the expected neutrino 
fluxes generated during propagation by the ultra-high en- 



ergy cosmic rays. We use recent results on the propaga- 
tion of UHE protons [go) to estimate the neutrino pro- 
duction. The aim is to explore the neutrino production 
with a photoproduction interaction model (SOPHIA [|T)) 
that fits well the experimentally measured multiparticle 
production data over a wide energy range. For this pur- 
pose we extend the calculation of proton propagation in 
the local universe by Stanev et al. [|o| to cosmological 
distances. We also study the importance of the cosmo- 
logical evolution of the sources of cosmic rays in different 
cosmological models. The aim of the present work is to 
study the level at which these ultra-high energy neutrinos 
are indeed guaranteed. 

Section O discusses the neutrino production from prop- 
agation of ultra- high energy protons in the local Universe. 
In Section III we obtain the neutrino spectra from homo- 



geneously distributed cosmic ray sources accounting for 
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the cosmological evolution of these sources. Section 
explores variations in this flux under the influence of dif- 
ferent assumptions concerning proton injection spectra, 
source evolution and distribution, and background cos- 
mology. Section [y| gives a brief overview of the event 
rates that could be expected in future large neutrino de- 
tectors. Discussion of the results and the conclusions 



from this research are given in Section VI 



II. NEUTRINO FLUXES FROM PROTON 
PROPAGATION IN THE LOCAL UNIVERSE 

We begin with our method for calculating neutrino pro- 
duction from proton propagation in the the nearby uni- 
verse (for a detailed discussion see [^(J ) . The calculation 
is carried out as a Monte Carlo simulation of individ- 
ual particle histories in the presence of the cosmic back- 
ground radiation, including energy loss processes such as 
photoproduction, e + e~ pair production, and adiabatic 
losses. The extension of this method to cosmological dis- 
tances is discussed Sec. 



III. 



An important new ingredient of the calculation is to 
use the event generator SOPHIA [^l| to simulate in de- 
tail the proton/neutron interactions with photons from 
the cosmic microwave background. This event generator 
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has several main differences from previously used codes: 

(1) the inclusion of direct pion production at the photo- 
production threshold. In this i-channel process the pho- 
ton is absorbed via a jtt + tt + vertex, and so only charged 
pions are produced. Although the cross section for this 
process is smaller than the dominant A + resonance, it 
yields a significant number of neutrinos, when folded with 
the steep proton injection spectrum. 

(2) explicit consideration of 10 different resonance pro- 
duction channels in the important energy region just 
above the particle production threshold. 

(3) QCD motivated multipion production at large center- 
of-mass energies. 

To evaluate the neutrino yields, protons are injected 
in narrow logarithmic bins and all products of their in- 
teractions are collected with the same energy binning. 
We use 10 bins per energy decade, and 10000 protons 
per bin weighted by an E~ 2 spectrum within each bin. 
The injection energy ranges from 10 19 to 10 23 eV. This 
gives the option to explore different injection power spec- 
tra and cutoff energies by rescaling the products of each 
energy bin. 
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FIG. 1. Neutrino production efficiency, summed over fla- 
vors, as a function of proton injection energy. The solid curve 
shows the ratio of the energy carried by neutrinos to that 
of electromagnetic particles due to photoproduction in fully 
developed cascades (200 Mpc), as simulated with SOPHIA 
[ ^l| . The dashed curve is the same ratio but for the A + reso- 
nance approximation which is frequently used in analytic cal- 
culations. The dotted curve shows the total neutrino energy 
relative to injected proton energy for SOPHIA. 

The results from the Monte Carlo are illustrated in 
Fig. |l|. Concerning the overall yield of neutrinos, the 
dominant feature is the turn on of the GZK process at 
E p w 5 x 10 19 eV. The ratio of yield in neutrino energy to 
yield in radiative energy depends primarily on the ratio 
of charged to neutral pion production. If all pion pro- 
duction occurred through the A + resonance this ratio 
would be approximately (| x \)/{\ + \ x \) — \, where 



for charged pions ~ | of the energy goes to neutrinos. 
At high energies, isospin "democracy" suggests that the 
ratio should tend to (| x |)/(| + \ x §) = 1. For low en- 
ergy protons, direct production of charged pions plays an 
important role, again increasing the neutrino yield above 
that expected from the A + resonance. 

The next step is to place the neutrino production 
model into an astrophysical setting. For the proton 
source spectra we use a power law with an exponential 
high-energy cutoff 

dN 

— cx E~ a xexp(-£/£ c ), (1) 

where a — 2 unless otherwise stated and E c = 10 215 
eV. During propagation, adiabatic energy losses for the 
protons are calculated assuming Hq = 75 km/s/Mpc. 
Similarly, neutrino energies are redshifted by a factor of 
(1 + z), where z is the redshift of the interaction site. 

The energy degradation of ultra-high energy protons 
in propagation in the microwave background is very fast. 
The minimum mean free path for photoproduction in- 
teractions is 3.8 Mpc at a proton energy of ^6xl0 20 eV. 
Protons with an energy of about 10 21 eV thus interact on 
the average twice or more during the first 10 Mpc of prop- 
agation and lose close to 50% of their injection energy. A 
significant fraction of the energy loss (about 40%) goes 
into neutrinos. The neutrino flux thus originates from 
the initial stages of proton propagation. 
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FIG. 2. Neutrino fluxes produced during the propagation 
of protons over 10, 20, 50, 100, and 200 Mpc (from bottom 
up) in a 1 nG random magnetic field. The heavy histogram 
shows the proton injection spectrum defined in Eq. ([j]). 

Fig. ^| shows the fluxes of electron and muon neutrinos 
after propagation over different distances up to a maxi- 
mum of 200 Mpc. About 60% of the final neutrino fluxes 
are generated in the first 50 Mpc and more than 80% in 
the first 100 Mpc. The contribution from the second half 
of the maximum propagation distance is small because 
the proton spectrum is deprived of > 10 20 eV particles 
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and photoproduction interactions are rare. It follows, 
therefore, that from the point of view of neutrino pro- 
duction a source at 200 Mpc produces a fully evolved 
spectrum. Accordingly, for the cosmological scenarios 

we scale the neutrino yields to 
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that follow in Section 
this result. 

There are two other noticeable features in the neu- 
trino spectra shown in Fig. |^. The muon neutrino spec- 
tra have a single peak at energies between 10 18 and 10 19 
eV. Electron neutrinos, however, exhibit a more compli- 
cated double peak structure. The first peak between 10 16 
and 10 17 eV is populated by v e from neutron decay. The 
neutron decay length equals the photoproduction inter- 
action length at about 4xl0 20 eV and neutrons of lower 
energy are more likely to decay than to interact. This 
leads to the formation of an additional peak in the elec- 
tron neutrino spectrum. The second peak, in a position 
similar to that of muon neutrinos, is populated mostly by 
v e from fi + decay with a small admixture of v e generated 
predominantly in neutron photoproduction interactions. 
The ratio of {v^ + v^) / {v e + v e ) in the second peak is 2, as 
expected, although integrated over the whole spectrum 
the ratio is closer to 1. 
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FIG. 3. Total number of neutrinos produced per injected 
proton of energy above 10 19 eV (lower curves) and 10 20 eV 
(upper curves). The proton energy was sampled from the 
spectrum ([!]). Solid lines show the sum of muon neutrinos 
and antineutrinos, the dashed lines are for electron neutrinos 
and antineutrinos. 

One also notices the slight shift of the peak of the dis- 
tribution to lower energy with the propagation distance. 
At longer propagation distance, lower energy protons suf- 
fer photoproduction interactions and generate lower en- 
ergy neutrinos. There is also a small effect from the adia- 
batic losses of all neutrinos, but this is hardly noticeable 
here because the maximum source redshift considered is 
only z w 0.05. 

Fig. |3| shows the total number of neutrinos produced 
per source proton as a function of the source distance. 
Because the protons lose most of their energy during 



propagation over the first 50 Mpc, one would naively ex- 
pect that the neutrino number does not change for source 
distances above 100 Mpc. The continuing rise of the neu- 
trino to proton number ratio is due to redshift effects. 
The minimum proton energy for photoproduction inter- 
actions decreases as (1 + z)^ 1 , which leads to an increase 
of the number of interacting protons. Even for relatively 
small redshifts involved (z — 0.1 for 400 Mpc) this leads 
to an increase of the generated number of neutrinos. 

Fig. |3| also underlines the importance of the maximum 
energy of the proton injection spectrum. In this calcula- 
tion we use an exponential cut-off with E c — 10 215 eV. 
Assuming an E c of 10 20 5 eV would not drastically de- 
crease the neutrino flux. Cutting off the proton injection 
spectrum at lower energy would, however, require very 
nearby sources for the extremely high energy showers de- 
tected by the Fly's Eye and AG AS A experiments. 



III. NEUTRINOS FROM PROTON 
PROPAGATION OVER COSMOLOGICAL 
DISTANCES 

In the following we will focus on the case of uniformly 
distributed sources with identical proton injection spec- 
tra. Although a homogeneous source distribution is dis- 
favored by the resulting source energy requirements and 
arrival proton spectra it serves here as a sim- 

ple generic model whose results can be easily rescaled to 
account for local density enhancements or even nearby 
point sources. 

The local neutrino flux of flavor i generated from the 
propagation of cosmic rays over cosmological distances 
can be written as an integral over redshift and the proton 
energy E p (s denotes "source") 



dEi 



C(z,E s p )Y(E s p: E Uil z)^dz. (2) 



Here, the neutrino yield function is 



Y(E°E Vi ,z) = E Ui 



dN p dE Vl 



Also, the source function per unit redshift is 
C(z,E' p ) =H(z)r,(z)C (E a p ), 



(3) 



(4) 



where H(z) parametrizes the cosmological source evo- 
lution, rj(z) describes the cosmological expansion, and 
Co(E p ) is a properly normalized version of the source 
spectrum in Eq. (Q). The metric element T){z) is defined 
as 



r](z) 



dt 

dz 



H (l 



Om(1 



+{l-Sl M -Sl A )(l + z) 



21-1/2 



(5) 
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which simplifies to 1/(Hq(1 + z) 5 / 2 ) for the Einstein - de 
Sitter universe (Hm = 1, = 0). 

The yield function Y is evaluated utilizing the Monte 
Carlo result for a 200 Mpc source and the scaling relation 



y{e;,e u , z ) = Y((i + z )e;, (i + z) 2 Pv,o). 



(0) 



In scaling E u one factor of (1 + z) arises from redshifting 
the neutrino energy from its observed value to its pro- 
duction value. Both E v and P* scale by (1 + z) to main- 
tain the same invariant reaction energies in the presence 
of a higher cosmic background temperature. Although 
it simplifies the numerical work considerably, utilizing 
the scaling relation introduces some approximations. For 
redshifts < 0.05 it overestimates the neutrino production 
as per Fig. ||. However the contribution to the total fluxes 
coming from z < 0.05 is very small (see Fig. [5] below). 
Another effect is that at high redshift the competition 
between neutron decay and neutron photoproduction is 
altered in favor of photoproduction, and so we make a 
modest overestimate of the v e flux around 10 16 eV. At 
high energies, the sum of u e and v e fluxes remains un- 
changed, but the flavor distribution may be altered. 
The source proton luminosity Cq is parametrized as 



C (Et) = P 




s dN p s 



psdNp 
p dE; ' 



(7) 



with dNp/dEp given by Eq. ([!]) and Po denoting the in- 
jection power per unit volume. 

The injection power of cosmic rays with energy above 
E m in = 10 19 eV can be roughly estimated using the 
local cosmic ray energy density p3]. The cosmic ray 
flux dN/(dEdQdAdt) at 10 19 eV is about 2.5xl0~ 28 
cm _2 s _1 ster~ 1 GeV _1 . Assuming that: 

1 ) all cosmic rays at that energy are extragalactic; 

2) 10 19 eV cosmic ray flux is as at injection; and 

3) the differential proton spectrum at injection is a power 
law with spectral index a — 2, 

one obtains a cosmic ray energy density 



Pe 



4tt 
c 



E 



dN 



dEdttdAdt 



dE 



(8) 



of 1.1 x 10 54 erg/Mpc 3 per decade of energy. To cal- 
culate the injection power required to maintain this en- 
ergy density one needs to make an assumption about the 
lifetime tcr of these cosmic rays. A conservative ap- 
proach would be to use a lifetime close to the Hubble 
time. Using tcr — 10 10 yrs gives a power of 1.1 xlO 44 
erg/Mpc 3 /yr per decade of energy. Of course, the total 
power for E >10 19 eV depends on the maximum energy 
at acceleration. 

The correct way of calculating the injection power for a 
model of cosmic ray source distribution and injection (ac- 
celeration) spectra is to propagate the accelerated spec- 
tra from the sources to us and fit the locally observed 
spectrum. We do not perform this procedure because it 



involves assumptions on the cosmic ray source distribu- 
tion and the structure and strength of the extragalactic 
magnetic fields which are beyond the scope of this paper. 
We use instead the cosmic ray injection power obtained 
in a similar, somewhat simplified way by Waxman pi] , 
who derived Po = 4.5 ± 1.5 x 10 44 erg/Mpc 3 /yr between 
10 19 and 10 21 eV for power law cosmic ray injection spec- 
tra with a between 1.8 and 2.7. We will use this value 
of Po for the energy spectrum of Eq. (|l|) integrated be- 
tween 10 19 and 10 22 eV. The higher P max approximately 
compensates for the factor of exp(— E/E c ) as compared 
to Waxman's result. 

Finally, we have to specify the cosmological evolution 
of the cosmic ray sources, TC(z). We use the parametriza- 
tion of [§31, i.e. 



H(z) = 



(1 + z) 3 

(1 + 1.9) 3 

(1 + 1.9) 3 exp{(2.7 



0/2.7} 



z < 1.9 

1.9 < z < 2.7 

z > 2.7 

(9) 



where n = 3 describes the source evolution up to mod- 
erate redshifts. We also briefly consider a stronger evo- 
lution model with n = 4 up to z=1.9 and flat at higher 
redshifts. 

Fig. |] shows electron and muon neutrino fluxes ob- 
tained with our nominal choice of astrophysical and cos- 
mological parameters, and carrying out the integration to 



a redshift of 



= 8. Integrating to infinity increases 



the neutrino fluxes by only about 5%. 

Fig. |^ also shows the limits on neutrino production in 



cosmic ray sources derived by Waxman & Bahcall [125 26 
(WB). As those calculations were carried out for the same 
source evolution model, similar spectra and the same in- 
jection power Po, they serve to compare the expecta- 
tions for 'source' versus 'propagation' neutrinos associ- 
ated with UHECR of astrophysical origin. Our propa- 
gation flux is slightly below the WB limit for the muon 
neutrino and antineutrino flux for energies between 10 18 
and 10 19 eV. The differences lie in the assumed neutrino 
yield per proton. For their limit, WB assume a maximal 
thin source, i.e. an energy equal to that of the injected 
proton is deposited into neutrinos, whereas for our cal- 
culation only a fraction goes into neutrinos, as shown in 
Fig. [|. At higher energies one can see the effect of the 
factor exp(—E/E c ) in our source spectrum. At lower en- 
ergies, the cosmic background radiation is devoid of high 
energy photons and so low energy protons do not pro- 
duce neutrinos. In contrast, the cosmic ray sources are 
assumed to have abundant higher energy photons and so 
the limit on source neutrinos continues to scale as P~ 2 
to low energies. 
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FIG. 4. Fluxes of electron neutrinos (dashed lines) and an- 
tineutrinos (dotted lines) generated in propagation of protons 
are shown in the upper panel. The lower panel shows the 
fluxes of muon neutrinos and antineutrinos. Solid lines show 
the sum of neutrinos and antineutrinos. The shaded band 
shows the Waxman & Bahcall |^,^6) limit for neutrino pro- 
duction in cosmic ray sources with the same injection power. 
The lower edge of the band is calculated without account for 
the cosmological evolution and the upper one with the evolu- 
tion of Eq. (0) . 

Fig. [5] is designed to show how the neutrino flux is 
built up from contributions at different redshifts. It is 
evident that the high and low ends of the neutrino spec- 
trum are sensitive to different epochs of the source evo- 
lution. First consider the protons that will contribute 
to neutrinos with energy 10 19 eV. At z — these pro- 
tons have an energy of a few times 10 20 eV, above the 
threshold for photoproduction. This energy will increase 
with the source redshift. As a result, the source contribu- 
tion EpdN/dEp for these neutrinos effectively decreases 
as (1 + z)~ x . To this we must add additional factors 
of rj(z)H.(z) ~ (1 + z) 5 for the source evolution in a 
Qm = 1 cosmology, and a factor of (1 + z) explicit in the 
(1 + z)d/d(l + z) plot. Together, the function plotted 
naively scales as (1 + z) 5 . This scaling stops at z = 1.9 
where Ti.(z) is assumed to flatten. For higher energy neu- 
trinos E v = 10 20 eV, the increasing proton energy runs 
into the exponential cutoff E c of our model injection spec- 



trum causing a further decrease with 1 + z. The result of 
these considerations is that the highest energy neutrinos 
are produced primarily by relatively young sources, and 
are sensitive to assumptions about the recent universe. 

For low energy neutrinos, say 10 16 eV the story is a 
bit more complicated. From kinematic arguments the 
prime production candidate for such neutrinos would be 
a proton of energy a few times 10 eV, but such protons 
are below the photoproduction threshold. Protons with 
higher energy can, of course, produce low energy neutri- 
nos, but due to the small phase space the production is 
suppressed by a factor of E v /E p . Now, as the source red- 
shift increases, E v at production also increases as 1 + z. 
At the same time, the minimum value for E p at produc- 
tion decreases due to the increasing cosmic microwave 
background temperature. Thus, phase space considera- 
tions of the neutrino production process yield a net factor 
of (1 + z) 2 . With the lowering of E p , the source spectrum 
factor yields an increase of 1 + z. Including r)(z)Ti(z) and 
the explicit 1 + z for the plot gives an overall dependence 
of (1 + z) 4 5 at low energies. This behavior continues un- 
til 

a) the source evolution model changes its z dependence, 
or 

b) the photoproduction threshold at z has dropped so 
that there is no phase space suppression for that neu- 
trino energy. At that point there is a transition to the 
high energy behavior outlined above. The net result of 
these considerations is that the low energy part of the 
spectrum is dominated by high redshift sources, and is 
sensitive to assumptions of a cosmological nature in our 
calculation. 
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FIG. 5. The curves, labeled by log 10 (-E„), show the contri- 
bution of different source distances to the neutrino flux as a 
function of redshift for our nominal n = 3 source evolution 
model given in Eq. (^). 

Finally, we comment on the energy where the neutrino 
flux peaks in Fig. ^. Given the turn on of photoproduc- 
tion (Fig. 0) and the kinematics of the A resonance, 
one might expect the peak to occur at around 10 19 eV. 
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Our Monte Carlo, however, yields more neutrinos with a 
softer spectrum than a A resonance model, so the peak 
from a low redshift source occurs at about 3 x 10 18 eV, 
as seen in Fig. 0. Moreover, as discussed just above, 
the peak of the cosmological spectrum is shifted by two 
factors of (1 + z) from the redshift which dominates the 
source contributions. For our H, r\ and £0 this occurs 
at (1 + z) = 2.9, and so the resultant neutrino spectrum 
peaks at around 3 x 10 17 eV as seen in Fig. |]. 

IV. VARIATIONS 

Many of the parameters associated with the calcula- 
tion of the neutrino fluxes shown in Fig. [| have rather 
large uncertainties. The power needed to maintain the 
flux of cosmic rays above 10 19 eV varies by about 30% 
for injection differential spectral indices between 1.8 and 
2.7 j24J]. The cosmological evolution of the source lumi- 
nosity evaluated from star formation regions ]27j ] could be 
somewhat stronger, as also indicated by the attempts to 
derive the cosmological evolution of GRB |^8| and their 
fluences [Q. We show the influence on the generated 
neutrino fluxes in Fig. ^ where we calculate the neutrino 
flux with the same injection power but a stronger cosmo- 
logical evolution - (1 + z) A up to z = 1.9 and constant 
thereafter. The stronger cosmological evolution increases 
the neutrino flux by a factor of 3 and generates a small 
shift of the maximum flux to lower energy. The integra- 
tion was carried again to redshift of 8. 



(crosses, energy cutoff at E = 10 215 eV). 

All three calculations show the peak of the neutrino 
spectrum at approximately the same energy of 2-3xl0 17 
eV. The spectrum of Yoshida & Teshima is somewhat 
narrower then the one obtained in this work, while the 
agreement with Protheroe & Johnson is very good. This 
latter work uses the cosmological evolution evolution 
model RLF2 @ with the 'fudge factor' of Rachen & 
Biermann ]3l{ ] rather than a simple redshift dependence. 

The neutrino flux calculated by Stecker et al. Jl^] (not 
shown) seems to be based on an injection power not much 
different from the normalization of Waxman. The spec- 
trum however peaks at higher energy. An error might 
have been made in accounting for the neutrino redshift 
(F.W. Stecker, private communication). 
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FIG. 6. The lower boundary of the shaded area corre- 
sponds to the neutrino flux shown in Fig. ^ with the source 
evolution n — 3 and the upper boundary is for n = 4 up to 
z — 1.9 and constant afterwards. Open diamonds show neu- 
trino fluxes calculated by Yoshida and Teshima Jl3| and the 
crosses are due to Protheroe and Johnson E3. 

The two sets of points in Fig. || represent calculations of 
Yoshida & Teshima (diamonds, source evolution with n 
= 4 and cutoff at z = 4) and of Protheroe & Johnson E3] 
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FIG. 7. Variation of neutrino flux resulting from different 
proton injection spectra. The dotted curve is for an injection 
spectrum of E~ 2 5 and the dashed curve is for E~ 3 . 

In Fig. Q we show the neutrino fluxes obtained with 
the source evolution of Eq. (||) and differential injection 
spectral indices of 2.5 and 3, keeping again the injection 
power and E c constant. The steeper injection spectra 
generate smaller neutrino fluxes at high energy, because 
of the much smaller number of protons above 10 20 eV 
that are mostly responsible for high energy neutrino pro- 
duction. By contrast, low energy neutrinos are predom- 
inately generated by protons injected at high redshift, 
where the threshold for photoproduction is decreased by 
a factor of 1 + z. Because the total injection power is 
kept constant, a steeper spectrum results in an increase 
in the flux of lower energy protons and, hence, low energy 
neutrinos. 

All results shown above are calculated for a homoge- 
neous source distribution. It has been suggested in the 
past, and recently in the context of a specific acceler- 
ation model in Ref. p^] that the observed cosmic ray 
spectrum can be best fit by a combination of a homoge- 
neous source distribution with an enhancement of local 
sources at distances less than 20 Mpc or with a single 



G 



source at a similar distance. 
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FIG. 8. Muon neutrino and antineutrino fluxes generated 
if one half of the local cosmic rays are generated by nearby 
sources at a distance of 20 Mpc. Note that fluxes are given 
for the full solid angle of 4-k. The dashed curve corresponds 
to one half of the muon neutrino flux of Fig. [l|, integrated over 
the full solid angle. The dotted line shows the production by 
locally generated cosmic rays by either a single source or by 
a local enhancement of the cosmic ray sources. 

Fig. H shows the neutrino fluxes generated under the 
assumption that the injection power of homogeneously 
distributed sources is Pn/2 and 50% of the UHECR at 
10 19 eV are generated by a single source at a distance of 
20 Mpc. This scenario will not predict a proton arrival 
spectrum similar to that of Waxman (24j] because the in- 
jection spectrum is kept, up to the normalization, the 
same for the single source at 20 Mpc and the homoge- 
neously distributed sources. It merely serves as a simple 
example of the changes to be expected in such a case. It 
is also a good approximation for a local density enhance- 
ment because the difference between the neutrino flux 
magnitudes for a single source and local enhancement 
scenarios is very small. 

The observational difference is substantial because, for 
a single source scenario, most of the high-energy neutri- 
nos would come from the direction of that source. Most 
of the neutrinos due to proton propagation from a lo- 
cal source would be generated in the first interaction of 
protons of energy above 10 20 eV (see Fig. ||). These high- 
energy protons do not scatter significantly in a random 
extragalactic field with an average strength of less than 
100 nG and the relativistic decay kinematics ensures that 
the neutrinos are emitted in the direction of the proton 
momentum. 

Finally we discuss the importance of the cosmological 
model. All calculations shown above are performed with 
the assumption of a flat, mass dominated universe (ftjf 
= 1). However, recent astrophysical observations agree 
better with models containing a cosmological constant 
A Q. From the behavior of 77 (see Eq. (g)) one may 



expect an increased contribution to the neutrino produc- 
tion from higher rcdshifts. 

Fig. [9] shows the difference in the expected neutrino 
fluxes for the Einstein - de Sitter Universe (solid line) 
and a model with CIm = 0.3 and fl\ = 0.7, as currently 
favored by measurements p4| , both with a source evo- 
lution proportional to (1 + z) 4 . We use the stronger 
source evolution and carry the integration out to a red- 
shift of 8 to emphasize the difference between the cos- 
mological models. For a flat universe the ratio of fluxes 
is always smaller than l/^l — Q\ (see Eq. (|])). The 
difference between the Einstein - de Sitter Universe and 
one with a non- vanishing cosmological constant does not, 
however, depend very strongly on the cosmological evo- 
lution model of the sources. The ratio between the two 
cosmologies is about 1.6 for n = 3 evolution and 1.7 for 
n = 4 evolution. 

We caution that in the absence of a model which con- 
sistently accounts for the effects of Ha on the source evo- 
lution function Tt(z), this increase should be regarded as 
an upper limit. Specifically, we have varied 77, but kept 
H constant, whereas it might be argued that keeping the 
product rfH constant would be a better approximation 
to the effects of a non-zero fl\ on the source evolution, 
in which case there would be no change in the neutrino 
flux. 
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FIG. 9. Muon neutrino and antineutrino fluxes from ho- 
mogeneous source distribution for the Einstein - de Sitter 
Universe (solid line) and for S7a = 0.7 (dashed line), using a 
cosmological evolution of (1 + z) 4 up to redshift of 1.9 and 
flat at higher redshifts. 

V. EVENT RATES FROM THESE NEUTRINOS 

Detection of neutrinos produced during the propaga- 
tion of UHECR protons is a challenge. The flux peaks 
above 10 eV where the neutrino nucleon cross-sections 
are of order 10 -31 cm 2 . Such values of <7„at are large 
enough to make the Earth opaque, but still require 100- 
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1000 km of water to ensure an interaction. As a result, 
there is no sensitivity to upward neutrinos and low ef- 
ficiency for downward neutrinos. Under these circum- 
stances the prime detector requirement is large mass, of 
order 100 km 3 (water), in order to guarantee a few events 
per year. With such large volumes, the typical event may 
be assumed to be a contained event where the visible en- 
ergy is dominated by the high energy shower associated 
with the neutrino interaction vertex. Special detector 
geometries may provide additional sensitivity to /i or t 
lcptons produced in charge current events. Although we 
phrase our discussion in terms of water/ice detectors, suf- 
ficient mass may also be achieved by monitoring large 
volumes of atmosphere. 

All types of neutrino interactions generate showers. In 
charged current (CC) interactions of electron neutrinos 
and antineutrinos the neutrino energy is completely re- 
leased in the form of shower particles. The hadronic 
shower carries a fraction y of the initial neutrino energy 
E v , while the electromagnetic shower carries the remain- 
ing (1 — y)E u . Although the electromagnetic shower will 
be stretched out by the LPM effect both showers 
are likely to be contained within the detector volume. In 
charged current interactions of muon neutrino and an- 
tineutrinos we only use yE v for the shower energy, as we 
do in the neutral current (NC) interactions of all neutrino 
types. We use the GRV98 ||3q] structure functions to cal- 
culate the neutrino cross sections at ultra-high energy. 
The results are similar to those in Gliick, Kretzer, and 
Reya |37) and to those calculated by Gandhi et al. ]3S|| 
and Kwiecinski, Martin & Stasto fl39|. 
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FIG. 10. Differential rates for showers initiated by charge 
current interactions of electron neutrinos (dashes) and an- 
tineutrinos (dots). The solid line gives the sum of those and 
the dash-dot lines gives the rate of the W~ resonance events. 
Absorption by the Earth is included so that at high energy 
the solid angle for detection is fi « 2n. 



by charged current interactions of electron neutrinos and 
antineutrinos for a cosmological source evolution with n 
= 3. The total shower rate is dominated by v e interac- 
tions around the peak energy of the neutrino flux. As the 
interaction cross section increases with energy, the signal 
is dominated by showers with energy above 10 8 GeV. At 
this energy there is no significant background, either from 
neutrinos or muons, produced by cosmic ray interactions 
in the atmosphere. The spike at E v = 6.3 x 10 6 GeV is 
the rate for W~ boson production by P e ~electron inter- 
actions p0| ], divided by 20. The differential rate in the 
vicinity of the resonance is very high but the relevant en- 
ergy range is small and the total rate is 0.03 per km 3 yr. 
Accounting for invisible W~ decay modes would reduce 
the rate for showers even further. 

Table [I| gives the shower rates per km 3 of water per 
yr for showers generated by different types of neutrino 
interactions and different flavors. These rates represent 
da v /dy folded with the flux of neutrinos reaching the de- 
tector and do not account for any experimental efficiency 
and detector biases. As such, they can only serve as an 
estimate of realistic event rates. The rates are low: the 
IceCube Q and ANTARES Q km 3 -size experiments 
may not expect to detect these neutrinos. In addition vi- 
olation of any of the following assumptions may further 
decrease expectations 

(i) the detected ultra-high energy cosmic rays are not 
produced by a single nearby source, 

(ii) the powerful sources of ultra-high energy cosmic rays 
are homogeneously distributed in the Universe, and 

(iii) our normalization of the power of the cosmic ray 
sources and the injection spectrum are correct. 

In case that only a fraction of the UHECR comes from ho- 
mogeneously distributed sources, these rates come down 
to the same fraction of the values presented in Table |. 
On the other hand, for the cosmological evolution model 
with n = 4 the rates will go up by a factor of about 3, and 
assuming a cosmological constant with Q\ — 0.7 could 
give another moderate increase by a factor of 1.6 - 1.7. 

TABLE I. Rates per km 3 water per year of showers above 
different energies generated by different types of neutrino in- 
teractions for Po = 4.5 x 10 44 erg/Mpc 3 /yr and a cosmological 
evolution with n = 3 for homogeneously distributed cosmic 
ray sources (see text). 
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all v, NC 


0.052 


0.046 


0.032 


0.008 


0.001 


v B , CC 


0.054 


0.051 


0.046 


0.024 


0.004 


+ u T , CC 


0.092 


0.080 


0.057 


0.014 


0.002 


total 


0.192 


0.177 


0.144 


0.046 


0.007 



Fig. nu shows the differential rate for showers initiated 
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The rates shown in Table Q do not depend on the as- 
sumption that the atmospheric neutrino anomaly could 
be explained by oscillations of into v T Q . On the very 
long astrophysical pathlength 50% of the muon neutrinos 
would be converted to tau neutrinos. Charge current in- 
teractions of v T would not, however, be different from the 
CC interactions and so their sum does not depend on 
oscillation parameters. 

Another interesting process is the detection of r's pro- 
duced in v T CC interactions in the material surrounding 
the detector. Estimates show that the tau decay inside 
the detector will generate a signal which is comparable 
in shower energy and rate to that of v e CC interactions. 
Whether or not r decays can be observed as separate 
events depends on the r energy, which affects both decay 
length and energy loss Q), and the detector volume. A 
proper study of event rates for throughgoing or stopping 
t and (i leptons depends on the detector geometry, loca- 
tion and surrounding material and is beyond the scope 
of this paper. 

VI. DISCUSSION AND CONCLUSIONS 

The biggest uncertainty in the magnitude of the neu- 
trino fluxes from proton propagation is related to the 
distribution of cosmic ray sources. All arguments about 
the normalization of the cosmic ray injection power and 
the cosmological evolution of the cosmic ray sources are 
based on the assumption that the detected UHECR are 
of astrophysical origin and are not accelerated at a single 
nearby source. If the latter were the case, and the nearby 
source were responsible for all particles above f 19 eV, we 
should have to restrict the source distance to less than 20 
Mpc |2(| . The local ultra- high energy cosmic ray density 
would then be much higher than the average cosmic ray 
density in the Universe. There will be, no doubts, other 
regions where ultra-high energy cosmic rays are acceler- 
ated and are over-abundant. The overall normalization 
of the cosmic ray power would then depend on the filling 
factor of such regions. If one estimates this filling factor 
from the volume of the walls of galactic concentration, 
such as the supergalactic plane, that describes the galac- 
tic distribution within a redshift less than 0.05 [45} ], to 
nearby voids it will certainly not exceed 10%. In such a 
case the neutrino production due to proton propagation 
would be minimal as shown in Fig. |^ with the dotted line. 

A correct estimate of the power in UHECR depends 
on the strength of the extragalactic magnetic field in our 
neighborhood. In case of an average strength of the tur- 
bulent field exceeding 1 nG, cosmic rays with an energy of 
about 10 19 eV and below (with a gyroradius of less than 
10 Mpc) will have diffusive propagation pattern, which 
will enhance their flux at Earth. On the other hand, reg- 
ular extragalactic fields may guide these particles along 
the walls of matter concentration. All effects related to 
the propagation of UHECR should be a subject of fur- 



ther investigation before we could give a more reliable 
estimate of the UHECR power. 

Recently several authors have argued that powerful 
astrophysical systems and potential cosmic ray sources, 
such as GRB, have a cosmological evolution stronger than 
(1 + z) 3 |^9| , ^6| . Rcf. presents a combined analysis 
of the far infrared luminosity as a tracer of the star for- 
mation rate |27|| as a function of redshift and the gamma 
ray burst fluence. The star formation rate is fit with 
an exponential rise of exp(2.9,z) to a redshift of 1.7 and 
the GRB fluence distribution suggests a slow decrease at 
higher z. The use of such strong cosmological evolution 
and high activity at redshifts higher than 2 would place 
the estimate of the neutrino flux halfway between the two 
cosmological evolution models shown in Fig. |9[ 

Another important factor is the injection spectrum of 
UHECR and the highest energy at acceleration. The 
qualitative picture of its influence is demonstrated in 
Fig. [|. Generally protons have to be accelerated to ener- 
gies above 10 20 eV to generate significant neutrino fluxes 
from their propagation. This threshold is reduced by 
(1 + z) for contributions from high redshift. In view of 
the observations of cosmic rays of energy significantly 
higher than 10 20 eV this is very likely if the sources of 
these particles are astrophysical objects. 

Assuming homogeneously distributed astrophysical 
sources we obtain a neutrino flux similar to the Waxman 
& Bahcall limit at neutrino energies above 10 18 eV. This 
means that, with the assumptions and restrictions dis- 
cussed above, one may expect similar fluxes of UHE neu- 
trinos produced in astrophysical sources and in UHECR 
propagation. 

Since the signature of these ultra-high energy neutri- 
nos are showers, different types of air shower detectors 
should also be able to observe the highest end of the 
neutrino spectrum. The effective volume of the Auger 
observatory for UHE neutrino interactions was estimated 
to 30 km 3 of water equivalent p7| . If this effective vol- 
ume was achieved for a shower energy of 10 19 eV, the 
Auger observatory would see about 0.3 events per year 
from the estimates shown in Table |. It is interesting 
to note that at an energy of 10 18 eV the r decay length 
(l T = 49E/ (10 18 eV) km) is of the order of the dimensions 
of the Auger observatory. 'Double bang' events, caused 
by t neutrinos as suggested by Learned & Pakvasa Q, 
could be detected if the sensitivity of the array was sig- 
nificant in this energy range. 

As the effective volumes required are extremely large, 
the proposed satellite air shower experiments EUSO and 
OWL §^,|5(| might be well suited for the observation of 
ultra-high energy neutrino fluxes; however, with the cur- 
rent advertised threshold of 5 x 10 19 eV most of the po- 
tential event rate would go undetected. Because of their 
large field of view, these detectors should in principle be 
able to observe 'double bang' events of energy above 10 19 
eV - if the threshold energy were lowered to that level. 

Showers generated by ultra-high energy neutrinos 
could also be observed by their radio emission |5l] , p2| . 







Prototype experiments are in operation |53| and sugges- 
tions have been made for full scale experiments | p4|]55| 
that would have an energy threshold of 10 18 eV and a 
full effective volume of 10 2 -10 4 km 3 . Such detectors could 
take advantage of the higher shower rate that corresponds 
to the maximum of the neutrino flux as shown in Fig. 

The potential detection of ultra-high energy neutrinos 
is a crucial experimental result that will help us distin- 
guish between an astrophysical (acceleration) and cosmo- 
logical (top-down) origin of UHECR. In top-down sce- 
narios the neutrino fluxes are primary, roughly equal to 
the gamma ray fluxes and at least an order of magni- 
tude above the ultra-high energy nucleon fluxes. In all 
astrophysical scenarios the neutrinos, due to cosmic ray 
interactions at their sources or in propagation, are sec- 
ondary and their flux is a fraction of the cosmic ray flux. 

Measuring neutrinos from CR propagation can also 
help to distinguish between protons and heavy nuclei, 
such as iron, as highest energy cosmic rays. The energy 
loss of heavy nuclei during propagation over cosmolog- 
ical distances is governed by photo-disintegration. The 
absorption of photons leads mainly to giant dipole reso- 
nance excitation of the nuclei and, with a high probabil- 
ity, subsequently to the emission of a single nucleon in the 
de-excitation process Q. Hence the neutrino spectrum 
is expected to be dominated by relatively low-energy neu- 
trinos {E v < 10 7 GeV) from the beta decay of neutrons 
and unstable nuclei. 

Finally it should be mentioned that the magnitude of 
the flux and the arrival direction of UHE neutrinos are a 
good indication of the cosmic ray source distribution in 
astrophysical scenarios. 
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